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Reduction of CCD Observations Made with the Fabry-Perot 
Scanning Interferometer. II. Additional Procedures 
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Abstract. We describe a software package used at the Special Astrophysical Observatory of the Russian Academy 
of Sciences to reduce and analyze the data obtained with the Fabry-Perot scanning interferometer. We already 
described most of the algorithms employed in our earlier Paper I (Moiseev. l2002|> . In this paper we focus on extra 
procedures required in the case of the use of a high-resolution Fabry-Perot interferometer: removal of ghosts and 
measurement of the velocity dispersion of ionized gas in galactic and extragalactic objects. 



1. Introduction 

Scanning Fabry-Perot interferometer (FPI) can be used 
to perform a detailed analysis of the structure and 
kinematics of galaxies, nebulae, and other extended ob- 
jects. Detailed descriptions of the main idea of this 
observational method and references to earlier papers 
can found in Moiseev (|2002p (hereafter Paper I) and 
Gordon et al. pOOOp . The result of observations has the 
form of a set of two-dimensional interferograms — the con- 
volution of the monochromatic image of the object with 
the transmission curve of the FPI at each step of scanning. 
After special reduction (phase-shift correction) these in- 
terferograms can be assembled into a "data cube" . In such 
a cube, each spatial element in the detector plane is asso- 
ciated with its individual spectrum. The spectral interval 
is usually not very wide and amounts to only 5-50 A, mak- 
ing it possible to study, e.g., the kinematics of ionized gas, 
based on the data for one or two emission lines. 

To reduce observational data obtained with the FPI, 
one needs appropriate software, which differs from the 
software employed to analyze the data obtained with slit 
spectrographs. Paper I gives a brief review of the systems 
commonly used to reduce such data. Of the recently pub- 
lished papers on the subject we point out the paper by 
Daigle et al. (|2006p who suggested a number of new pro- 
cedures: adaptive spatial binning of data cubes prior to 
constructing the velocity fields, algorithms for detecting 
emission lines and subtraction of night-sky lines. 

At the Special Astrophysical Observatory of the 
Russian Academy of Sciences (SAO RAS) the scanning 
FPI is a part of the CCD-based SCORPIO multimode 
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focal reducer (Afanasiev & Moiseev, I2005p . Unlike two- 
dimensional photon counters, which made a return re- 
cently and are now used in combination with the FPI 
(Gach et al.. [2002p . "slower" CCDs require a different kind 
of algorithms to reduce the data obtained. This concerns, 
first and foremost, photometric correction of channels and 
night-sky spectrum subtraction. We addressed this issue 
in detail in our earlier Paper I, where we described the 
algorithms used in the software to reduce SCORPIO ob- 
servations performed in the scanning FPI mode. 

According to ASPIE0 database, a total of about 180 
FPI data cubes have been obtained with the 6-m telescope 
of SAO RAS in 2000-2007, with at least 30 publications 
based on these data. The experience in the reduction of 
this observational data allowed us to improve our software 
package and develop a number of new useful procedures 
for data analysis. In this paper we focus on describing 
these procedures. In section [2] we address the problem of 
removing the ghost light, which arises in high-resolution 
FPIs, in section [3] we consider the problem of velocity dis- 
persion measurements and in section |4] briefly describe 
the software package currently employed to reduce obser- 
vations made with the 6-m telescope using the scanning 
FPI incorporated into SCORPIO instrument. 

2. Ghosts subtraction 

SCORPIO is now equipped with two scanning interferom- 
eters (hereafter referred to as IFP235 and IFP501) with 
a gap between the plates corresponding to 235 and 501 
orders of interference, respectively, at the wavelength of 
A6562.8A. See Paper I for detailed parameters of both 
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interferometers. Because of the FPI location in the out- 
put pupil of the system between the collimator and the 
camera of the focal reducer ghost lights appear on the 
image. Ghost lights are due to backlight reflecting inside 
the interferometer plates and between these plates and the 
nearest lenses of the focal reducer, surfaces of the narrow- 
band filter, etc. This is a well known problem, and var- 
ious ghost "families" were described by Bland-Hawthorn 
(|1995p . and a somewhat more detailed description can be 
found in the paper by Jones et al. (fSOOS]) . Putting IFP235 
into the beam produces ghosts so that the main image 
and these ghosts are aligned symmetrically with respect 
to the optical axis of the FPI. The intensity of these ghost 
reflections (which we refer to as "diametral ghosts" , "D" , 
in accordance with the ghost classification suggested by 
Bland-Hawthorn (|1995p and Jones et al.. 12002]) amounts 
to about 10% of the intensity of the main image. In such 
cases it is usually recommended to tilt the FPI with re- 
spect to the optical axis of the system so as to move the 
ghost outside the detector field. However, certain design 
features of SCORPIO focal reducer prevent the use of this 
option. Therefore during observations the image of the ob- 
ject should be placed off the optical axis, and care must 
be taken to ensure that ghosts from field stars do not fall 
on the object of interest. All these problems decrease the 
working field of view of SCORPIO focal reducer by a fac- 
tor of two. 

In observations made with IFP501 another — 
"exponential" ("E", according to the terminology of 
Bland-Hawthorn PMS]) and Jones et al, [^UU^ — familv 
of ghost images appears. These ghost images form as 
a result of back reflections inside the plates of the FPI 
and pose a more serious problem than diametral ghosts, 
because the nearest exponential ghost is located within 
only 16" of the object image, and its intensity is rather 
high (about 12% of the brightness of the object). Such 
intense ghost images are due to the degradation of the 
antireflective coating of the outer surfaces of the IFP501 
plates. Unfortunately, the cost of plate replacement (or 
recoating) in this case is close to that of assembling a 
new FPI, which is rather high. We nevertheless were 
tempted to use IFP501 before a new interferometer was 
acquired, because the former performs well in terms of 
other characteristics — first and foremost, because of its 
R « 8300 spectral resolution at the Hq, line wavelength, 
which is rather high for extragalactic astronomy. Such 
an instrument proved to be very popular among the ob- 
servers at the 6-m telescope for tasks that are impossible 
to perform with IFP235 because of its lower spectral 
resolution. These were primarily observations of nebulae 
with relatively small (within 100-200 km/s) range of 
radial velocities: Herbig-Haro objects, star-forming 
regions in the Milky Way and other galaxies, etc. 

In a number of cases ghost traces can be removed at 
the stage of data reduction by applying the appropriate 
algorithms. Figure[l] (whose idea we borrowed from Bland- 
Hawthorn (|1995p and Jones et al. , '2002!) clearly illustrates 
the complex pattern of ghost images. Namely, the image 



of the galaxy (G) produces closely-spaced ghost images 
EqG and EiG whose brightness is equal to 4% and 12% 
of the object brightness, respectively. Now, EiG produces 
E2G, which, in turn, produces E3G — and the brightness of 
these ghosts is equal to 1.4% and 0.15% of the brightness 
of G, i.e., the brightness of ghosts decreases exponentially 
as it must be in the case of multiple reflections. These 
ghosts appear at the front plate of the FPI before light 
passes through the interferometer, and therefore are re- 
ferred to as monochromatic — their relative intensity in the 
data cube does not change with wavelength A. At the same 
time, ghost EqG appears at the back plate of the FPI after 
interference has taken place. The image of such a ghost on 
the interferogram is located in a domain corresponding to 
other wavelengths. Hence the data cube should exhibit an 
offset in A between the spectra of G and EqG, and it must 
depend on the position of the image on the detector. 

In Fig. [1] we also indicate diametral ghosts, i.e., ghosts 
located symmetrically with respect to the optical axis of 
the interferometer. The image of galaxy G and the corre- 
sponding ghost DC are located symmetrically with respect 
to the optical axis, ghost EiG produces diametral ghost 
DEiG, and ghost E2G, respectively, produces DE2G, etc. 
Ghost DG also generates a family of secondary exponen- 
tial ghosts: EiDG can be seen in the figure. The latter is 
superposed by diametral ghost DEqG. Diametral ghosts 
have rather complex structure, as is immediately appar- 
ent in the case of ghost DS produced by the brightest 
foreground star. Namely, this ghost consists of two images 
of the star — a normal image (with the intensity equal to 
^ 5% of that of star S) and a strongly defocused image, 
which has the shape of a ring, whose total intensity is 
equal to about 10% of intensity S. 

Note that it is by no means always possible to find 
so many ghost images as in Fig. [T] The figure illustrates 
a limit case corresponding to observations of a relatively 
bright object, whereas usually only the brightest ghosts 
(DG, EqG, and EiG) can be seen. We describe the pat- 
tern of ghost images in such detail, because we believe 
that our analysis will help the users better understand 
the data obtained and distinguish real emission features 
from ghosts. 

If we denote the ideal data cube (i.e., the cube without 
ghosts) as Ireai, then the wavelength-calibrated observed 
cube has the following form: 

Iobs{x,y,\) ^ Ireaiix,y, \) + Ighost{x,y, X), (1) 

where the family of exponential ghosts can be represented 
in the form: 

Ighost{x, y, A) w foPi;^}^i{x - Axq, y - Ayo, z)+ 

+ J2i=T filreal {x - Axi , y - Ayi , A) . 

Relative brightness fi of ghosts decreases exponen- 
tially with the number i of reflections and therefore sum- 
mation is performed not to infinity, but up to imax = 3—5, 
because the contribution of higher-order ghosts is usually 
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Fig. 1. Examples of ghost images produced in the IFP501 interferometer. Ha observations were performed on January 
30/31, 2007 at the request of C. Muhoz-Tuhon. Here "G" and "S" denote the II Zw 70 galaxy and a bright foreground 
star, respectively. The sum of all channels in the data cube is shown. Diametral and exponential ghost images are 
denoted by "D" and "E" , respectively. The cross indicates the center o the field of view. 



insignificant. For the same reason, we neglect the contri- 
bution of secondary ghosts, such as non-monochromatic 
satellites of the ghosts of the monochromatic family. The 
first term in ([2]) describes the wavelength-shifted ghost 
EqG. Here P denotes "phase correction", i.e., the transfor- 
mation from spectra expressed in terms of the FPI chan- 
nel numbers z to spectra expressed in terms of the wave- 
lengths: 

y, A) = PI{x, y, z) = I{x, y, k{z + p{x, y))), 

where p{x, y) is the so-called phase map (see Section 2 in 
Paper I for details), and k{z) = kiz + fcp is a linear func- 
tion whose coefficients are determined by the particular 
FPI employed. Correspondingly, /^J^; is a result of the re- 
verse transformation from wavelength scale to the scale of 
interferometer channels: 

Our aim is to infer from lobs the best approximation 
to ghost-cleaned cube Ireai ■ We use the following iterative 



procedure for this. We substitute lots for Ireai in formula 
(|2| to obtain the first approximation to the ghost model 
Igfiost(-^reai) ■ In this casc, the first approximation to the 
ghost-cleaned cube has, according to formula the fol- 
lowing form: I^^^; = lobs - IIhosv We now substitute ll^^i 
into formula ^ to derive the following approximation to 
ghost model Igi^ost^ ^^c, up to Igf^ost- Because of the small 
relative brightness of ghost images {fi <C 1), only a few 
iterations (n =3-4) are needed to construct a close-to- 
reality ghost-image model. 

The following relations are true for the family of 
monochromatic ghost images (i > 1): 

Axi = iAxi, Ayi = lAj/i, fi = f, 

and hence model ^ has only six free parameters, which 
are chosen in a way to minimize the differences between 
^J/host ^1^*^ ^obs ■ The following table gives the mean values 
of these parameters, which vary little from night to night. 

Our experience in the reduction of the data obtained 
with IFP501 shows that the model based on the above 
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Fig. 2. Subtraction of ghost images. Here we show the sum of aU channels of the data cube for the II Zw 70 galaxy: 
(a) initial image, (b) model of brightness distribution in ghost images, and (c) the result of subtracting the model 
from the initial cube. 



Table 1. Parameters of ghost-image model in Ha line 
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algorithm usually describes real ghost images rather well, 
with no appreciable bias. The only exception are the data 
obtained under unstable atmospheric conditions, when 
counts in individual channels had to be averaged in order 
to compensate for appreciable (more than 10-20%) varia- 
tions of the FWHM of stellar images (for a description of 
photometric correction see Paper I. 

Figure [2] shows how the ghosts around the image of 
the dwarf galaxy II Zw 70 are subtracted. Other exam- 
ples of observational data processed using the above de- 
ghostification procedure can be found, e.g., in Lozinskaya 
et al. (2 003.2006 ) - a study of star-forming regions in the 
dwarf galaxies IC1613 and VII Zw403, Martmez-Delgado 
et al. (|2007p - kinematics of ionized gas in blue compact 
galaxies, Movsessian et al. (|2007p - study of outflows from 
young stellar objects in the HL/HX Tau region. 

Yet another problem arises in observations of objects 
with strong surface-brightness gradients, where an intense 
ghost of, e.g., the galaxy nucleus, projects onto regions of 
much lower brightness. The level of Poisson noise is de- 
termined by the combined intensity of the ghost and the 
base. Therefore after the subtraction of the ghost model 
a situation may arise where the useful signal in the region 
considered is comparable to the amplitude of photon noise. 
To avoid loss of spectral information from low-brightness 
regions in observations of such objects, we recommend to 
divide the planned exposure into two and perform obser- 
vations successively with two different orientations of the 
instrument's field of view turned by about 90° in posi- 
tion angle. The ghosts obtained in the two data cubes 
should then project onto different regions of the object. 
After primary reduction, the corresponding ghost model 
is subtracted from each data cube and the regions are 



masked where the signal-to-noise ratio decreased strongly 
after ghost removal. The two cubes are then combined 
into one and spectra of masked regions in each cube are 
substituted by the corresponding "good" spectra from the 
other data set. One of the best examples of application of 
the above algorithm is the reduction of observations of 
the nearby dwarf galaxy IC 10 reported by Lozinskaya et 
al. (|2008p . Here the image of emission shells of ionized 
gas occupies more than half of the entire field of view of 
SCORPIO instrument. Figure [3] illustrates the sequence of 
operations in the process of subtracting ghosts from the 
images of these galaxy. 

3. Velocity dispersion measurement 

A number of observational programs performed with the 
FPI (e.g., mapping the ionized-gas velocity dispersion in 
galaxies) require accurate estimates of the halfwidths of 
emission lines. These estimates must take into account 
the broadening due to the instrumental profile. It is com- 
mon practice to use the following formula (which assumes 
that both the instrumental profile of the spectrograph and 
the initial — unbroadened — line profile can be described by 
Gaussian functions with the dispersions equal to areai and 
cTins, respectively): 



Cobs 



- . U2 



(3) 



Here aobs denotes the dispersion of the Gaussian used to 
describe the profile of the line observed at the output of 
the instrument. Hereafter by measuring the velocity dis- 
persion we mean estimating (Treat from the observed spec- 
tra. One should keep in mind that broadening of lines in 
the spectrum of the observed object may be caused not 
only by velocity dispersion agas (the measure of chaotic 
motions along the line of sight), but also by a number of 
other factors. Thus, according to Rozas et al. (|2000p . the 
following formula can be written for the integrated spec- 
trum of HII regions: 



real 



'gas 



'N 



(4) 
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Fig. 3. Sum of channels in the data cubes for the IC 10 galaxy. Ha-line observations were made on September 8/9, 
2005 at the request of T. A. Lozinskaya. The figure shows a 72" x 72" fragment of the image, (a) Initial image in the 
first cube; (b) the image after ghost subtraction, the arrow indicates regions with unsatisfactory quality of subtraction; 
(c) the image in the cube obtained by turning the field of view by 90°; (d) subtraction of ghosts from the second cube 
(the arrows indicate low-quality regions), and (e) the combinations of the two data cubes. 
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where cr^ « 3 km/s and a^^ « 9.1 km/s correspond to 
the natural width of the emission hne and its thermal 
broadening at 10'' K, respectively. 

Because of its simple and self-explanatory form, for- 
mula ^ is widely used to analyze spectroscopic data. It is 
often generalized by substituting FWHM for cr. However, 
one must keep in mind that exact equality in ([3]) is 
achieved only for Gaussian functions, because a convolu- 
tion of two Gaussians is also a Gaussian. The assumption 
about the Gaussian form of the instrumental profile is usu- 
ally true for slit spectrographs. However, the instrumental 
profile of the FPI, which is given by the Airy function, has 
wide wings and is best approximated by a Lorentz pro- 
file (see, e.g., Bland-Hawthorn (|1995p and Paper I) rather 
than by a Gaussian. Therefore if the initial profile of the 
emission line is a Gaussian with the dispersion determined 
by formula ^ , then the observed profile is the convolution 
of the Gaussian and Lorentz profiles and is hence given by 
the Voigt function: 

v{x y) = 1 y r ^~"'dx 

where 

A — Aq Wins 

a — —= , y = —= . 

\J 2a real V^Ureal 

Here Aq is the central wavelength and Wms denotes the 
halfwidth of the instrumental (Lorentz) profile of the FPI 
as determined from the spectrum of the lines of the cal- 
ibration lamp. We then approximate the observed profile 
by function ([5]) to obtain the required estimate areai- It 
is evident from Fig. [3] that compared to the Gaussian the 
Voigt profile fits much better the line wings in the galaxy 
spectra observed with IFP501. It is clear from general con- 
siderations that formula is a more correct tool for es- 
timating the velocity dispersion than formula ([31) ■ In the 
latter case we have to use Gaussian approximation for the 
profiles of the lines that deviate systematically from the 
adopted approximation (Fig.U^), and this may introduce 
an additional error in the estimated u^eai- However, this 
approach is highly popular in velocity-dispersion measure- 
ments for extragalactic HII regions (see, e.g., Martinez- 
Delgado, [20071 Muhoz-Tuiion, [l995l Rozas et al. 120001) . 
We believe that this method owes its popularity not only 
to the less complex appearance of formula Q compared 
to that of formula ©J but also to the fact that Gaussian 
approximation of spectral-line profiles is incorporated into 
virtually all packages of astronomical data reduction. The 
use of the Voigt profile for analyzing extragalactic spec- 
tra is less common (see, e.g.. I1994j) . despite the fact that 
integration in ([5]) poses no problem for modern computers. 

We estimated the errors of measurement of the kine- 
matical parameters (radial velocity and velocity disper- 
sion) for both approaches considered. We smoothed the 
instrumental profile of the FPI by a Gaussian with the 
dispersion equal to Gin and then added noise to the result- 
ing spectrum and estimated the velocity dispersion using 
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Fig. 4. Example of the Ha-line spectral profile in the 
II Zw 70 galaxy based on observations with IFP501 (dots). 
The solid line shows Gaussian (a) and Voigt (b) fits. 



both methods. The difference between the output {a out) 
and input {a in) velocity dispersions allows us to estimate 
the error aerr of the velocity dispersion measuring. We ob- 
tained a total of 1000 independent measurements for each 
fixed signal-to- noise ratio [S/N). We similarly estimated 
the error of measured radial velocity. We performed our 
computations for the instrumental profiles with the width 
near the Ha line equal to Wins — 35 and Wins = 115 km/s 
for IFP501 and IFP235, respectively. Figures [5](a, b) and 
[ni(a, b) show the results of computations — radial-velocity 
and velocity-dispersion errors as functions of the signal 
level. As expected, the error of radial-velocity measure- 
ments for symmetric lines does not depend on the algo- 
rithm employed. At the S/N — 30 it is equal to 2.5 and 8 
km/s for interferometers IFP501 and IFP235, respectively. 

The situation is quite different for the error of the mea- 
sured velocity dispersion. As expected, the error of mea- 
surements based on Voigt-profile fits smoothly decreases 
with decreasing signal level and amounts to only several 
km/s for S/N > 20. These errors are due to noise in the 
measured spectra and they contain no systematic compo- 
nent. Contrariwise, if measurements are based on relation 
(PI) then systematic but not random component of mea- 
surement error starts to dominate at S/N > 10. More 
specifically, velocity dispersion is overestimated, as it is ev- 
ident from Figs. [BJ; and [5};. Velocity-dispersion estimates 
inferred from Gaussian fits to the profiles exceed the ac- 
tual values by 7-8 and 20-25 km/s for IFP501 and IFP235, 
respectively. Such a systematic error is unimportant for es- 
timating widths of lines with areai > 100 km/s. However, 
velocity dispersion may be overestimated by up to 100% in 
the case areai = 10 — 20 km/s. This may be of critical im- 
portance, e.g., in the studies of ionized gas in star-forming 
regions, when it is necessary to identify expanding shells 
or regions where the velocities of chaotic motions exceed 
the speed of sound in the interstellar medium (Marti'nez- 
Delgado et al. [20071 Relafio & Beckman. [20051) . In these 
cases velocity dispersion estimates inferred by fitting the 
Voigt profile are to be used. Or, if for some reasons the 
authors prefer Gaussian approximation, it is necessary to 
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Fig. 5. Simulation of errors of measurement of kinematic parameters (at the la level). The computations were made 
for IFP501. The diamond signs and filled circles show the Gaussian and Voigt-profile approximations, respectively, (a) 
Dependence of the error of radial-velocity measurements on the signal-to-noise ratio. The average velocity dispersion 
was equal to 50 km/s; (b) error of measured velocity dispersion as a function of signal-to-noise ratio under the same 
conditions; (c) comparison of the initial and measured velocity dispersion (for S/N = 20). 




Fig. 6. Same as Fig. [5l but for IFP235. Figures (a) and (b) are for the case where an average velocity dispersion is 
equal to 200 km/s. 



estimate systematic errors like we did it above and correct 
correspondingly the estimates based on formula ([3]). 

Note that Relano & Beckman (W05) proposed an al- 
ternative method to account for the instrumental profile 
of the FPI using the reconstruction (deconvolution) tech- 
nique. This procedure is to be used for multicomponent 
lines. However, the method can be applied only to spectra 
with sufficiently high signal-to- noise ratios. 



4. Brief description of the software package 

To reduce observations made with the FPI on SCORPIO, 
we wrote IFPWID program package with multi- 
window user-friendly interface in IDL 6.X language 
(Fig. [7]). The codes are publicly avail- 
able at|http://www.sao.ru/hq/moisav/soft. Below in 
this section we describe the principal sequence of data- 
reduction steps to be made when using these programs. 
We do not describe in detail most of the algorithms 
employed — wavelength-scale calibration, photometric cor- 
rection, night-sky line subtraction, etc. (see our sufficiently 



detailed Paper I, and only briefly list the initial data for 
the reduction of a set of CCD frames: 



— OBJECT — inteferograms of the object studied. 

— NEON — images of interference rings from the emis- 
sion line selected by the narrow-band filter from the 
spectrum of the He-Ne-Ar calibration lamp. This cal- 
ibration is usually performed before and after the ob- 
serving night. 

— FLAT — interferograms of uniform "flat-field" illumi- 
nation produced by the continuous-spectrum lamp and 
obtained with the same narrow-band filter as was used 
for observations of the object. 

— TEST — images of the rings from the line-spectrum 
lamp in some FPI channels obtained along with ob- 
ject integrations at the same position of the telescope. 
These images are used to monitor the accuracy of scan- 
ning and to control the offset of the center of the rings 
due to the instrument flexures. 
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Fig. 7. IFPWID interface: the main menu (on the left) and the menu of operations with caHbrated cubes (on the 
right). 



4.1. Assembling the Data Cubes 

The first stage of reduction usually requires no special 
settings. As a rule, it suffices to indicate only the num- 
ber of the observational night, file-name templates, and 
the required frame format (some observations with the 
FPI involve the readout of only a fragment of the CCD). 
The program extracts all the remaining necessary informa- 
tion (parameters of the FPI, numbers of spectral channels, 
etc.) from the descriptors of the corresponding FITS files, 
which are filled automatically in the process of observa- 
tions. Therefore the user only has to check the required 
reduction steps in the menu (see the upper part of the 
menu in Fig. [7]) and press the "run" button. Below we 
briefly list these steps (their names used in the reduction 
program are italicized): 

— Search for and averaging of the bias-current frames 
taken in the required mode of CCD readout (MeanBias 
creation). The resulting super bias frame is then sub- 
tracted from all object integration and calibration 
frames. The CCD employed contains virtually no "hot" 
pixels, dark current is small and therefore can be ne- 
glected for exposures about several minutes. 

— Creation of data cubes from individual frames {Cubes 
creation). Superbias is subtracted and bad columns are 
masked. Further operations are the reduction of three- 
dimensional cubes (object interferograms and calibra- 



tion data assembled in order of channels). Cubes are 
stored in the standard FITS format (NAXIS=3). 

— Removal of cosmic-ray hits from the cubes of 
calibrating-lamp integration with line (NEON clean) 
and continuous (FLAT clean) spectra. Simple a-fHter 
is used here: the counts in the spectra that deviate from 
the mean by more than preset threshold value are sub- 
stituted by the half-sum of the neighboring channels. 

— Construction of the phase-shift map based on the re- 
sults of Lorentz-profile fits to the lines in the NEON 
cube (Phase map). A correction is applied, where nec- 
essary, to allow for nonuniform scanning of the cali- 
brating cube (see Paper I). 

— Testing the accuracy of wavelength scale (Phase test). 
The NEON cube is corrected for the phase shift and 
then the position of the emission line of the calibration 
lamp is measured for each pixel. 

— Computation of corrections to the wavelength scale 
constructed from the calibration cube (Test rings). 
The offsets (along both coordinates in the CCD plane 
and along the wavelength coordinate) of TEST frames 
are computed relative to the NEON cube (to correct 
for the instrument flexures and to monitor the scan- 
ning accuracy of the FPI). 

— Correction (if needed) of relative variations of FLAT 
lamp brightness during scanning of the calibration 
cube (Flat norm). 
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— Measurement of the offsets between the transmission 
maxima of the narrow-band filter measured for the 
FLAT cube and for the spectra of selected stars in the 
field of the object {Flat Z-shift). This procedure is usu- 
ally needed only for observations made with IFP235, 
where the half-width of the narrow-band filter is ap- 
preciably smaller than the wavelength interval between 
the neighboring orders of interference. The offset is 
usually close to zero, however, it may reach 2— 3.5 A for 
some filters. 

— Division of the cube by the flat-field (to correct for 
the transmission curve of the narrow-band filter) and 
cosmic-ray hit removal (OBJ flat and clean). 

4.2. Subtraction of Night-Sky Lines and Conversion to 
the Wavelength Scale 

Then comes the turn of the reduction procedures that re- 
quire the user's intervention more often than during the 
stage of assembling the data cubes. The procedures corre- 
spond to the buttons of the main interface (the left-hand 
part of Fig. [7^), which are italicized below. Reduction is 
usually performed in the following order: 

— Correction for the background gradient {2D gradient). 
In case of observations of objects in some narrow-band 
filters the background brightness distribution in the 
object cube exhibits appreciable gradient even after 
the division by the FLAT, especially in the presence 
of additional light pollution due to the Moon. This 
residual gradient must be due to nonuniform illumina- 
tion from the FLAT lamp combined with specifics of 
the interference coatings of particular filters. The user 
can correct this effect by setting the parameters of the 
two-dimensional brightness distribution to which the 
object cube is normalized. 

— Creation of the mask for subtracting the spectrum of 
the night sky {Mask OH). The mask is based on the 
image of the sum of the channels of the object cube. 
Regions with the brightness below the given threshold 
are considered as "background" . The resulting mask 
can be edited if necessary. 

— Setting the parameters for the sky-background sub- 
traction {SKY parameters). We described our tech- 
nique of sky-background subtraction in the Paper L 
Its main idea is to average sky background emission in 
each channel of the OBJ cube over the azimuthal angle 
within narrow rings centered on the optical axis of the 
FPI with taking into account the mask constructed at 
the previous stage of reduction. The average bright- 
ness profile is then subtracted from the interferogram 
of the object. The user may choose various modes of 
such subtraction, vary the width of the rings where 
averaging is performed, fix the center of the rings or 
set options for an automatic search of the ring center 
in each channel. When needed, averaging can be per- 
formed within individual sectors (such a procedure can 



be useful for correcting for variations of the instrumen- 
tal contour over the field of view), etc. 

— Subtraction of the night-sky spectrum in accordance 
with the parameters given above {Sky remove). 

— Photometry of stars in each channel of the object cube 
( Channels preparation) . The relative offsets of the im- 
age centers, variations of FWHM and integral flux are 
measured for stars from a precompiled list and the av- 
eraged dependences of these parameters on the number 
of the channel in the cube are constructed. 

— Photometric correction of the object cube based on 
the results of photometry of stars in each individual 
channel {Photo/Shift corrections). Corrective channel 
offsets are applied and account is taken of variations 
of atmospheric transparency (the channel counts are 
multiplied by the computed coefficients) and seeing 
(convolution with two-dimensional Gaussians). 

— Conversion of the object cube to the wavelength scale 
{Linearization) . 

— Subtraction of ghosts using the procedure described in 
Section [2] ( G/ioste remove). 

4.3. Processing of Calibrated Cubes 

The reduction steps described above produce a 
wavelength-calibrated data cube with maximum ac- 
count taken of all instrumental effects. In this form it can 
be analyzed using various software tools depending on the 
user preferences or on the task to be accomplished. The 
data can also be converted to the format adopted in the 
popular ADHOCH reduction system. Below we describe 
the operating sequence to be performed with our software 
tools in order to map emission-line radial-velocity and 
velocity-dispersion distributions . 

Figure [7] (on the right) shows the menu of operations 
with calibrated cubes. The following procedures (their 
names are italicized) are performed: 

— Rotation of the cube to the "correct" orientation of the 
images (with North at the top and East on the left), 
because observations can be performed at any arbi- 
trary (or specially selected) position angle {Rotation 
to NE). The turn angle is computed from the data of 
the descriptors of the FITS-file header to within 0.1°. 
In the cases where better accuracy is required astro- 
metric reduction is to be performed using field stars. 

— Correcting the radial velocities for the motion of the 
Earth about the Sun {V -heliocentric). The necessary 
information is extracted from the FITS-file header. 

— If necessary, a cube fragment is cut containing the ob- 
ject studied {Sub-cube). In this form reduced FPI cubes 
are usually stored in our ASPID database [I] 

— Smoothing the cubes using a one- and two- 
dimensional Gaussians of given width along the spec- 

^ ADHOC software package was developed by 
J. Boulesteix (Observatoire de Marseille) and is available 
at http : / /www . oamp . f r/ adhoc/ adhocw . htm 
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Fig. 8. Interface of the program of the analysis of emission-hne profiles. 



tral [Smoothing along Z) and spatial {Smoothing along 
X) coordinates, respectively. 



We construct the velocity fields, velocity-dispersion 
maps, monochromatic images in the emission line and in 
the continuum, using GAUS program whose main inter- 
face is shown in Fig. [5] This program can be used to ana- 
lyze individual profiles in selected pixels of the data cube 
by fitting the emission-line profiles to Voigt and Gauss 
functions. The integral in ([5]) is computed via standard 
VOIGT function of the IDL language. The main task per- 
formed by the program is an automatic identification of 
emission lines in the data cube and their approximation 
by various functions for the given parameters. This is im- 
portant because the number of pixels that contain useful 
signal in our data cubes may amount to several hundred 
thousands, thereby preventing individual approach to each 
spectrum. Our procedure yields two-dimensional maps of 
profile parameters (Doppler velocities, line fluxes, etc.). 



5. Conclusions 

In this paper we briefly describe the software currently 
used to reduce observations made with the scanning FPI 
operated as a part of SCORPIO instrument. We believe 
the software complex described here to fully meet the re- 
quirements imposed by the research tasks performed at 
the 6-m telescope of the SAO RAS using the observa- 
tional technique considered. In our opinion, further evo- 
lution of the data reduction software package should be 
associated with project of the development of SCORPIO- 
2 new-generation multimode focal reducer at the SAO 
RAS. We hope, first and foremost, that the new illu- 
minator of the calibration beam would make it possi- 
ble to abandon flat-field correction procedures (nonuni- 
form residual background, FLAT cube wavelength shift 
observed in some narrow-band filters etc.). Higher au- 
tomation level of SCORPIO-2 (compared to the current 
version of the instrument) will allow further unification of 
the process of acquisition data with the FPI and manage 
without "manual" setting of parameters in a number of 
procedures. The latter primarily concerns subtraction of 
the night-sky spectrum. In this case it will suffice to simply 
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develop a pipeline for the reduction of the data obtained 
with the scanning FPI. 
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